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Abstract 

We study dynamics near the threshold for blowup in the focusing nonlinear Klein-Gordon equa- 
tion uu — u xx + u — \u\ 2a u = on the line. Using mixed numerical and analytical methods we find 
that solutions starting from even initial data, fine-tuned to the threshold, are trapped by the static 
solution S for intermediate times. The details of trapping are shown to depend on the power a, 
namely, we observe fast convergence to S for a > 1, slow convergence for a = 1, and very slow 
(if any) convergence for < a < 1. Our findings are complementary with respect to the recent 
rigorous analysis of the same problem (for a > 2) by Krieger, Nakanishi, and Schlag [3]. 
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INTRODUCTION 



In this paper we consider the nonlinear Klein-Gordon equation (NLKG) with a focusing 
power nonlinearity in one spatial dimension 



for real-valued u(t,x) and smooth compactly supported (or exponentially localized) even 
initial data. This equation arises in a variety of physical applications, for instance as a 
model of self-focusing waves in nonlinear optics. It is well known that for some initial data 
solutions of Eq.(l) are globally regular in time and for some initial data they blow up in 
finite time. In particular, the global regularity for small data follows from standard energy 
estimates [1] , while blowup for negative energy data follows from a concavity argument [2] . 
In addition, not all globally regular solutions decay to zero, as is evident from the fact that 
there exists the (unstable) static solution 



The dichotomy of global regularity vs blowup raises two obvious questions: (i) what 
determines a borderline between these two behaviours and (ii) what is the evolution of 
critical initial data which lie on the borderline. These questions have been recently studied 
by Krieger, Nakanishi, and Schlag [3]. They proved that there exist a codimension-one 
center-stable manifold associated with S which locally separates the regions of blowup and 
dispersion, and solutions lying on this manifold are trapped by S (i.e., remain in its small 
neighborhood) for all future times. Moreover, they proved that for a > 2 trapping actually 
implies asymptotic convergence to S, however the rate of convergence was not determined. 

The aim of this paper is to obtain a precise quantitative description of the asymptotic 
dynamics on the center-stable manifold (decay rates and spatial profiles). To this end, 
we proceed in two steps using a combination of analytical and numerical methods. In 
step one, we derive analytically the asymptotic behaviour of linearized perturbations about 
S. In step two, we employ the numerical bisection method to fine-tune initial data to the 
borderline between blowup and global regularity. Such specially prepared solutions approach 
S, hang around it for some time (whose span depends on the accuracy of fine-tuning), and 
eventually are ejected out along the one-dimensional unstable manifold of S. We find that 



u tt - u xx + u - \u\ a u = 0, 



a > , 



(1) 




(2) 
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this behaviour is very well approximated by the linearized dynamics obtained in step one. 
The main outcome of our study is that the dynamics on the center-stable manifold of S 
exhibits an interesting dependence on the power a which is a consequence of the change 
of spectral properties of the linearized operator around S as a varies. More precisely, we 
distinguish three different scenarios of convergence to S as t — > oo: (a) fast dispersive decay 
\\u(t, x) — S(x)\\ oo ~ t ^' 2 sin(t) for ol ^ 1, (b) slow dispersive decay H^x^t, x) — *5'(^')||oo ^ 
t _1 ' 2 sin(t) for a = 1 due to the presence of the zero energy resonance, and (c) very slow (if 
any) decay for < a < 1 due to the presence of oscillatory modes. 

The paper is organized as follows. In section 2 we recall some well-known facts about the 
spectrum of the linearized operator around S. These facts are used in section 3 to describe 
the linearized dynamics near S. In section 4 we present numerical simulations of the critical 
dynamics for Eq.(l) for three representative powers, a = and |, and confront them 
with the results of section 3. Finally, in section 5 we mention several open problems for 
future work. 



SPECTRUM OF THE LINEARIZED OPERATOR 



In order to make the paper self-contained and to fix notation, in this section we recall the 
well-known facts about the spectrum of the linearized operator around the static solution S. 

Substituting u(t,x) = S(x) + f(t,x) into Eq.(l) and linearizing, we obtain the linear 
evolution equation for small perturbations around S 

fu-f xx + f + V(x)f = 0, V{x) = - {2a + l) J 2 a + l) . (3) 

cosh ax 

After separation of variables, f(t, x) = e st v(x), we get the eigenvalue problem 

Lv = -\ 2 v, L = -JP. + V (x), (4) 

where A 2 = s 2 + 1. The potential V(x) is the exactly solvable Poschl- Teller potential. Since 
we are interested only in even eigenf unctions, we restrict the domain of the operator L 
to square-integrable functions on the positive half-line with the Neumann boundary condi- 
tion v'(0) = 0. Then L is self-adjoint with the continuous spectrum [0, oo) and a discrete 
spectrum depending on a. The generalized eigenf unctions have the form 

v {a \X,x) = (cosh axf +1/a) 2 F 1 ( 1 + 1 + i; - sinh 2 a^ . (5) 

\ 2a 2a 2 J 
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Using the connection formula for hypergeometric functions (for a non-integer a — b) 
2 F l {a, b, c; z) = ~ °| (1 - z)-° 2 Fi ( a, c - 6, 1 + a - 6; 



r(6)r( c -a) v ' V i-* 

we rewrite (5) in the form (for A 7^ 0) 

v {a) (x,\) = A a (\)v { «\\,x) + A a (-\)v { +\\,x) , (7) 

where 

r(i)r(A) 

^ (A) = r(i + ^)rH-^)' (8) 

1±A 1 1 =F A A 2 \ , n , 
1 ± — ; cosh ax . (9) 



4 a) (A,x) = {coshax)^ a 2 Fx(l + 



{a \\,x) ~2 ±x /°e TXx , (10) 



2a 2 2a a 

For x —7- 00 we have 

^± 

so u_ '(A, x) and u+ (A, a;) represent ingoing and outgoing waves, respectively. Eigenvalues 
are given by the positive roots of the equation A a (\) = 0: 

A n (a) = a + 1 - 2na, n = 0,1, . . . N , (11) 

where N is the largest integer less than (a + l)/2a. Note that due to the dispersion relation 
A 2 = s 2 + 1, eigenvalues with A > 1 correspond to unstable modes which grow exponentially 
in time as e^ x2 ~ lt , while eigenvalues lying in the interval < A < 1 correspond to oscillatory 



(neutral) modes behaving as sin(\/l — A 2 1). Thus, for each a there is exactly one unstable 
eigenmode (with the eigenvalue Ao(a) = a + 1), while for a < 1 there are, in addition, N 
oscillatory eigenmodes. Denoting by v^f 1 (x) an eigenfunction associated with the eigenvalue 
A n (a) we get from (7), (10), and (11) 

vi a \x) := v^\\ n ,x) = A a (-\ n )v ( °:\\ n ,x) (12) 

= (coshax)~^ 1+1 ^ a ' )+2n 2F1 ( — — n H , — n, 2 — 2n H ; cosh~ 2 ax 

\2 a a 

The hypergeometric function above is a polynomial in cosh _2 ax of order n. Factorizing this 
polynomial we get 

n 

v V)(x) = (coshax)" (1+1/Q) Y[ (1 - c nk sinh 2 ax) , (13) 

k=i 



with certain positive coefficients c n k- For example, the first two eigenfunctions are 

v { Q a \x) = (coshax)- {1+1/a) , X ( ^ = a + 1, (14) 

v[ a \x) = (coshax)- (1+1/Q) (l - ^ sinh 2 a^ , \{ a) = 1 - a . (15) 

An important role in our analysis will be played by the generalized eigenfunction at the 
endpoint of the continuous spectrum A = 

v ( a \0,x) = (coshaxf +1/a) 2 Fi ( 1 + 1 + l ; -sinh 2 ax J . (16) 

\ 2a 2a 2 / 

Using the formula (see 15.8.9 in [4]) 

jp f \ r(c)(l - z)~ a ^ (a) fc (c - a) k k 

2Fl(fl ' a ' ^ Z) = T(a)T(c-a) g (jfeip (1 " z) 

x (ln(l -z) + 2ifj(k + 1) - ^(a + k) - ip(c - a + k)) , 

we obtain the asymptotic behaviour for x — > oo 

v (a) (0, *) = — -j^— -r- x + (9(1) . (18) 

If (a + l)/2a is a positive integer (i.e., a = 1, 1/3, 1/5, ...), then the coefficient of the leading 
order term in (18) vanishes and therefore v^ a \0,x) is everywhere bounded, which means 
that there is a resonance at the endpoint of the continuous spectrum (this is also seen as 
the zero "eigenvalue" in (11)). 

Remark 1. Note that negative zeros of the coefficient ^4 Q (A) correspond to antibound states. 
For a ^ 1/k (k G N) they are given by two infinite series 

A" (a) = a + 1 - 2na and A+ (a) = -{2a + 1 + 2ma) , (19) 

where n runs over all integers greater than (a + l)/2a and m runs over all non-negative 
integers. For a = 1/k there are no antibound states because of the cancelation of infinities 
in the gamma functions in the numerator and the denominator of A a (X). 

Remark 2. In two important cases of quadratic and cubic nonlinearities the generalized 
eigenfunctions j)W(A,i) can be expressed by elementary functions. For a = 1/2 we have 

(A 2 - iyi)(A,x) = cosh(Ax) (a 2 + H _ (20) 

- tanh (!) f 3 A 2 + \ > (21) 

A \2J I 4 8cosh 2 f J ' V ; 



and for a = 1 we have 



(A 2 - l)v (1) (X,x) = cosh(Ax) (A 2 + 2 -=- J - 3Atanh(x) sinh(Ax) . (22) 

V cosh x J 

LINEARIZED DYNAMICS 

We will show in the next section that one can prepare special initial data for which the 
solution approaches the static solution S and remains close to it for some time. It is natural 
to expect that during this transient phase of evolution the dynamics can be approximated 
by the linearization around S which leads us to the study of Eq.(3). For t > the solutions 
of this linear equation are given by 

f(t,x) = j d t G(t,x,y)f(0,y)dy + j G(t,x,y)d t f(0,y)dy, (23) 

where G(t, x, y) is the retarded Green function. In order to determine the asymptotic be- 
haviour of Green's function for late times we shall use the Laplace transform in time 

POO 

g(s,x,y) := / e- st G(t,x,y)dt. (24) 
Jo 



which is well-defined and analytic in s for Res > y a(2 + a). Following the standard 

procedure, we construct g(s, x, y) from two linearly independent solutions 1 of Eq.(4), v (A, x) 

and v + (X,x), satisfying the appropriate boundary conditions: v'(\,0) = and t> + (A,x) ~ 
e -\x f or x _j, ^ 

/ % ~ / \ \ 1 \ v(\,x)v+(\,y) for x < y 

g{s : x,y) =g(\,x,y) := — — t (25) 

W \ A ) \ v{\,y)v + (\,x) for y < x 

where W(X) is the Wronskian of solutions v(\,x) and v + (X,x). Taking advantage of the 
fact that the Wronskian does not depend on x and computing it at infinity using (7) and 
(10) we get W(X) = -2XA(X). 

The function g(X,x,y) is analytic in A for Re (A) > Ao and has a meromorphic continua- 
tion to the whole complex plane of A with a countable number of poles located at the zeros 
of the Wronskian. The function g(s,x,y) inherits all these poles through the dispersion 



relation A = \/ s 2 + 1 and, in addition, has two branch points at s = ±i. We take the branch 



1 Hereafter, to avoid notational clutter we drop the superscript (a) on the eigenfunctions. 
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cuts emanating from the points s = ±z and going horizontally to the left (Im (s) = ±1 and 
Re (s) — oo). This choice is convenient because only the poles s n corresponding to the 
eigenvalues y/s^ + l = A n > given by (11) lie on the first Riemann sheet while the poles 
corresponding to the antibound states A^ < given by (19) lie on the second sheet 2 . 
Next, having g(s,x,y), we compute G(t,x,y) via the inverse Laplace transform 

i pc+ioo 

G{t,x,y) = — / g(s,x,y)e st ds, (26) 

where c is taken to the right of all the singularities, which in our case means c > a(2 + a). 
In order to obtain the late-time asymptotics t ^> \x — y\ of G(t, x, y), we deform the contour 
of integration to the left half-plane as shown in Fig. 1 and use the Cauchy residue theorem. 




FIG. 1: The deformation of the contour of integration. 

For late times t ^> \x— y\ the integral over the semi-circle at infinity vanishes so G = Gr+Gc, 
where Gr is the contribution from the sum over residua at the poles of g(s, x, y) lying on the 
first sheet and Gc is the contribution from the integrals surrounding the two branch cuts. 
We have 

N N 

G R = y2 Res[fi'(s,x,y)e si ] = a e Sot v (x)v (y) + V]a n sin(w n t + S n )v n (x)v n (y) , (27) 

n=0 n=l 

where s = \J a(a + 2) and oo n := a/1 — A 2 with A n given in (11). 

2 With this choice of cuts, for a pole A > 1 there are two poles s = ±a/A 2 — 1 on the first sheet. As A 
decreases and crosses the value +1 these poles meet at s = and continue as a pair s = ±iyl — A 2 . 
For A = they reach the branch points ±i. Decreasing A below zero causes the two poles go back to 
s = ±ivl — A 2 and further to s = ± \/A 2 — 1 but now on the second sheet. 
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The integrals along the branch cuts, because of the exponential damping e st for Re (s) < 0, 
are dominated by contributions from the endpoints s = ±i. The behaviour of g(s, x, y) near 
the endpoints follows from the Laurent expansion of g(X,x,y) at A = 

g(X, x, y) = "-1^1 + Co(x , y) + Cl(x , y)X + 0{\ 2 ) , (28) 

where c_i is nonzero if and only if there a resonance at zero energy. Let A ± g(s) be the jump 
of g(s,x,y) across the upper/lower cut. From (28) we get A ± g(s) ~ (s =p i) 13 near s = ±i, 
where (3 = —1/2 if there is a zero energy resonance and otherwise /3 = 1/2. This yields the 
late-time asymptotic behaviour 

1 r t^ 1 / 2 sin(t) c_i(x, y) resonance at zero, 

G c = ~ g{8,x,y)^da^i (29) 
Jcuts \t 3 / 2 sin(t) c\(x, y) no resonance at zero. 

It is evident from the expression (25) that c_i(x, y) ~ v(0,x)v(0,y). Somewhat surprisingly, 
if zero is not a resonance, the same is true for Ci(x, y). To see this, let us differentiate Eq.(4) 
with respect to A and take the limit A — > 0. We get L(d\g)(0, x, y) — — limA_>o 2A^(A, x, y) = 
0, hence d\g(0,x,y) is a zero mode of L. Since by construction d\g(0,x,y) satisfies the 
Neumann boundary condition at zero and is symmetric in x, y, it follows that c±(x, y) ~ 
v(0,x)v(0,y), as claimed. 



NUMERICAL RESULTS 

In this section we present numerical simulations of the critical behaviour for Eq.(l) on 
the half-line x > with the Neumann boundary condition ut(t, 0) = and exponentially 
localized initial data. The numerical technique is standard. We use the method of lines with 
a fourth-order Runge-Kutta time integration and fourth-order spatial finite differences. All 
simulations were performed in the quadruple (128-bit) arithmetic precision. 

Before presenting the results let us recall that Eq.(l) has the conserved energy 

E(u, ut) = I f°° (u 2 t +u 2 x + u 2 - \u\ 2a+2 \ dx . (30) 

2 Jo V a + 1 / 

The static solution S is the critical point of the static energy functional E (u) := E(u,Q). 
Another useful functional is 

j poo 

K(u):=±E (J3u)\ fh=1 = / (u 2 x + u 2 -\u\ 2a+2 )dx. (31) 
dp J 
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Since SE (S) = 0, it follows that K(S) = 0. The significance of the functional K was first 
pointed out by Payne and Sattinger [5] who showed that for solutions with E(u, Ut) < Eq(S) 
one has dichotomy: global regularity for K(u) > and blowup for K(u) < 0. Note that 
the boundaries of the Payne-Sattinger sets /C+ = {(u,u t )\E(u,u t ) < Eq(S), K(u) > 0} and 
JC- = {(u,u t )\E(u,Ut) < E (S), K(u) < 0} are separated (except for a single point (S, 0)), 
hence the borderline between global regularity and blowup necessarily lies above E (S). 
The results shown below correspond to the following one-parameter family of initial data 3 

u(0,x) = (a + l)^exp , u t (0,x) = 0, (32) 

Note that these data lie entirely above the energy E (S). We find that solutions with small 
a are globally regular while solutions with large a blow up in finite time. Using bisection 
we find that there is a single a* which separates these two regions. 

We remark in passing that by repeating the bisection for many different families of initial 
data (depending on one, two, or even three parameters) one can probe the shape of the 
center-stable manifold in the phase space. Recently, very interesting studies in this direction 
were reported in [6], however we do not pursue them here since our goal is different. We 
want to get the precise quantitative description of the critical evolution. For this purpose, 
as emphasized above, the choice of an interpolating family of initial data is irrelevant. 

At cr* the curve of initial data (32) intersects the center-stable manifold of S. Of course, 
numerically it is impossible to start the evolution exactly at this intersection point, so 
hereafter we show only near-critical solutions (with a* determined to about 30 decimal 
places). Such solutions, by continuity, approach S and stay in its neighborhood for some 
time until they are ejected out along the one-dimensional unstable manifold of S. On the 
basis of linearized analysis from section 3 we expect that when the near-critical solution 
stays close to 5* it should be well approximated by the expression 

N 

u(t, x) ~ S(x) + A e Sot v (x) + ^A n sm(oo n t + 5 n ) v n (x) + C r 1 sin(t) u(0, x) , (33) 

n=l 

where so = ^/a(a + 2), oo n = ^1 — (A n ) 2 , and 7 = 1/2 (if there is a zero energy resonance) 
or 7 = 3/2 (otherwise). For critical data the coefficient Ao(a*) vanishes, however in practice 

3 We stress that the near-critical dynamics is universal in the sense that it does not depend on a particular 
choice of a family of initial data which interpolates between basins of global regularity and blowup. 
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there is always a small admixture of the unstable mode which means that the length of time 
during which we can observe the near-critical dynamics is limited by a characteristic time of 
growth of the instability ~ — ^ In \ Sa\, where 5a is the distance from the threshold achieved 
by bisection. 

The results of our numerical simulations are depicted in Fig. 2,3,4 for three powers a = 
1/2, 1, and 3/2 (corresponding to quadratic, cubic, and quartic nonlinearities, respectively). 
They illustrate three qualitatively different scenarios of critical dynamics. In each case we 
show two plots: (a) time evolution at a fixed point in space (for convenience we take x = 0) 
for a pair of marginally sub- and super-critical initial data (32) with a = a* ± 10~ 30 , and 
(b) snapshots of spatial profiles at several times before ejection. Below, in order to facilitate 
the comparison of numerics with the predictions of linear theory, we give the explicit form 
of the expression (33) specialized to each case. 

• a = 3/2: The linear approximation (33) takes the form 

u{t, x) ~ S(x) + A t v (x) + C t~ 3/2 sin(t + 8)v(0 7 x) , (34) 



where 



S(x) 



v (x) 



v(0, x) 



, 3x 
cosh — 
2 



3.r 



( cosh 
V 2 



( cosh 
V 2 



„ 8 

3oA 3 



4 4 1 . , 2 3x 
-,-,-: — smh — 
3 3 2 2 



(35) 
(36) 
(37) 



1.5 



1.45 



1.35 



1.25 





FIG. 2: a = 3/2. (a) Fast dispersive decay to S, (b) Convergence of the spatial profile to the zero 
energy eigenfunction (37). The critical behaviour for all a > 1 is qualitatively the same. 
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a — 1: The linear approximation (33) takes the form 

u(t, x) ~ S{x) + A e^voix) + Ct~ 1/2 sin(t + 6) v(0, x) , (38) 

where 

S(x) = \/2cosh _1 (x) , (39) 

vq(x) = cosh~ 2 (x) , (40) 

v(0,x) = 3cosh" 2 (x) - 2. (41) 



1.5 
1.48 
1.46 
1.44 
1.42 

1.4 
1.38 
1.36 
1.34 



subcritical 
supercritical 
sqrt(2) + Csin(t+8)/t" 




1 

0.5 



j -0.5 

-1.5 

-2 
-2.5 




2 3 4 
X 



FIG. 3: a = 1. (a) Slow dispersive decay to 5, (b) Convergence of the spatial profile to the zero 
energy (resonance) eigenfunction (41). 



a = 1/2: The linear approximation (33) takes the form 



(t, x) ~ S(x) + A ^o(x) + A 1 sin(^t + 5x) v^x) + C t~ 3/2 sin(t + 5) y (0, x) , (42) 



u 



where 



S(x) 
v (x) 
vi(x) 
v(0, x) 



of X 

2( COSh 2 

(cosh |) 3 

x ^ ~3 
cosh — 
2 

15 / x 
— cosh — 
4 V 2 



1-4 sinh' 



.r 



-2 11 



hi tanh ( — 

4 V2 



3 15 / , i\- 
cosh — 

4 8V 2 



(43) 
(44) 
(45) 
(46) 
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FIG. 4: a = 1/2. (a) Long-lived oscillation around S, (b) Convergence of the spatial profile to the 
oscillatory eigenmode (45). 

The above results confirm that the near-critical dynamics is well described by the lin- 
earized approximation (33). If there are no oscillatory modes (which happens for a > 1), 
then we observe asymptotic convergence to S in the form of an algebraic oscillatory tail 
with frequency 1 and the decay rate 7 = 3/2 for a > 1 and 7 = 1/2 for a — 1. The spatial 
profile of the tail is shown to converge to the zero energy eigenfunction. We emphasize that 
the oscillatory tail should not be confused with the so called quasinormal ringing which is 
governed by a different mechanism [7]. 

For a < 1 the late-time dynamics is dominated by the oscillatory modes. At the linear 
level the critical solution is asymptotically periodic (if there is only one oscillatory mode, 
which happens for l/3<a<l)or quasiperiodic (if there are two or more oscillatory modes, 
which happens for < a < 1/3). We expect that the nonlinearity will induce a slow decay 
of the oscillatory modes via the resonant transfer of energy to radiative modes, however the 
calculation of this nonlinear decay is not an easy task, both numerically and analytically. 
Numerically, because even for very accurate bisection with 5a « 1CT 30 , the "lifetime" of 
near-critical dynamics is much too short to perform a reliable measurement of decay of 
amplitudes of oscillatory modes. Analytically, because the linear dispersive decay seems to 
be too weak to allow the calculation of radiation damping along the lines of [8]. Note that 
for a = 1/2 we have 2u\ > 1, hence the frequency 2uj\ lies in the continuous spectrum, which 
indicates that the nonlinear decay should occur already at the second order perturbation 
level. In Fig. 5 we verify that the deviation from the linearized approximation does in fact 
oscillate with frequency 2oj\. 
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0.001 
0.0005 


-0.0005 
-0.001 



-0.0015 

20 25 30 35 40 45 

t 

FIG. 5: a = 1/2. The deviation from the linearized approximation (obtained by subtracting the 
linear fit from the numerical solution) is shown to oscillate with the frequency 2oj\. 

FINAL REMARKS 

We hope that the explorations of the asymptotic behaviour of critical solutions described 
in this paper will inspire other researchers to carry on the investigation and prove our 
conjectures. We conclude with several remarks. 

1. In this paper we considered only even initial data. This restriction considerably facil- 
itates the analysis since it eliminates translation and boost symmetries and prevents 
the solution S to move. For general data these symmetries must be accounted for by 
the modulation analysis. We believe that this is only a technical difficulty and our re- 
sults for even data hold true in the general case provided that the assertions about the 
conditional asymptotic stability of S are replaced by the conditional orbital asymp- 
totic stability (with suitably red-shifted frequencies of oscillatory tails and modes) (cf. 
[9, 10] for the related problem in three space dimenions). 

2. We take this opportunity to point out that the asymptotic behaviour of solutions of 
Eq.(l) which decay to zero is not quite well understood. For a > 1 the globally 
regular small data solutions are asymptotically free for t — > oo, however for a < 1 
the asymptotic effect of the nonlinearity cannot be neglected and the solutions do not 
scatter to free solutions [11]. In the latter case one has to analyze the long range 
modulation of solutions, which (to our knowledge) has been accomplished only for 
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a = 1 [12]. Without the smallness assumption, as far as we know, the dynamics of 
globally regular solutions in one space dimension is an uncharted territory (in higher 
dimensions the dispersion is stronger and the problem becomes easier, cf. [13, 14]). For 
instance, it would be interesting to explore the post-ejection dynamics of marginally 
subcritical solutions and see how they enter into the small data regime. 

3. As we saw above, perturbing the NLW equation by the mass term significantly affects 
the dynamics of globally regular solutions. In contrast, the asymptotics of blowup 
is structurally stable under this perturbation. Both for the NLW and the NLKG 
equations the blow-up is governed by the ordinary differential equation u tt — \u\ 2a u = 0, 
hence u ~ (T — t)~ 1 / a as t /*■ T. This was first proved for the NLW equation by 
Merle and Zaag [15] and recently generalized to a class of perturbed NLW equations 
(including NLKG) by Hamza and Zaag [16] (in one space dimension these results hold 
for all a > 0). 

Acknowledgments: The authors are grateful to the Erwin Schrodinger Institute in Vienna, 
where part of this work was done in February 2010 during the workshop "Quantitative 
Studies of Nonlinear Wave Phenomena". PB acknowledges an inspiring discussion with 
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